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' The infra-red properties of three-dimensional abelian lattice gauge theory are known to be governed 

by a neutral plasma of magnetic monopole excitations. We address the fate of these monopoles in 
the presence of light dynamical fermions, using a lattice formulation of compact QED3 with Nf — 4 
fermion flavors supplemented by a four-fermi contact term permitting numerical Monte Carlo simu- 
lations in the chiral limit. Our data hint at a restoration of chiral symmetry above a critical value of 
the (inverse) coupling /3. By performing simulations in a sector of non-vanishing magnetic charge, 
we are able to study the response of the theory to an external magnetic test charge. Our results sug- 
' gest that the monopole plasma persists even once chiral symmetry is restored, and hence survives the 

continuum limit. 

o ' 

O ' 1 Introduction 



finement in QCD.^ His particular sense of beauty in physics has brought to the field the 
idea that confinement can be understood in terms of a symmetry . The role of magnetic 



Adriano's outstanding interest in recent years has been an understanding of color con- 

<D , . .. .... ^ 

^ monopoles can then be exposed by constructing a disorder parameter that tests the mag- 

netic properties of the vacuum.- Many of the concepts can be exported to other theories 
where dual excitations are conjectured to play a role in determining the phase structure of 
the system; compact abelian lattice gauge theory in 3 dimensions is one such theory. This 
model, which we will refer to as "quenched compact QED3", has been understood for 
many years as a field theory example where interesting IR effects such as linear charge 
confinement and generation of a mass gap arise as an effect of monopoles™^! (which 
in 3d are not particles but instantons). As the continuum limit is approached the system 
remains in the confining phase found at strong coupling, with the vacuum thought of as 
a dilute charge-neutral plasma of monopoles m and antimonopoles m whose effect is to 
Debye-screen the Coulomb potential of a test charge. The photon mass is proportional 
to exp(— C/e 2 ), e being the electric charge appearing in the relation between gauge po- 
tential and link variables U^ x = exp(iaeA Atx ), where a is the lattice spacing and C(a) 
related to the free energy of an isolated Dirac monopole. Although this scale is 0(a~ 1 ) 
implying that the monopole has a core extending over a few lattice spacings where UV ar- 
tifacts are significant, following Polyakov's original treatment^ the model can be viewed 
as an effective description of the Higgs phase of a 3d Georgi-Glashow model in which 
semi-classical solutions exist for m and rh, and hence a continuum limit identified as 
= l/e 2 a — >■ 00. 



More recently people have begun to ask if this situation persists if light electrically- 
charged fermions are present, in which case the model is now "compact QED3". A simple 
way of seeing why this might be expected to have a dramatic effect is to observe that as 
a result of the Dirac quantisation condition the combination eg, where g is the magnetic 
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charge of the monopole, is a renormalisation group invariant^ 1 implying g^ > g. Hence, 
virtual electron anti-electron pairs anti-screen the mm interaction. More detailed cal- 
culations^ * suggest that V m fh(x) is modified from the 3d Coulomb form 1/x to lnx, 
implying that a deconfined phase with monopoles only present within tightly-bound mm 
molecules, and therefore unable to influence the IR physics, may exist for sufficiently 
large (3. Other workers, however, have argued that interaction among magnetic dipole 
pairs restores the Coulomb potential at large x, and that confinement continues to survive 

the continuum limitlr^ 

A related question is whether there exists a chiral symmetry restoring transition at 
some finite C (like confinement, spontaneous chiral symmetry breaking leading to a mass 

gap for fermion excitations is known to happen at strong coupling — ). In non-compact 
QED3 in which the Wilson action is replaced by a simple F 2 term which is a non- 
periodic function of the plaquette, so that monopoles are suppressed and play no role in 
the continuum limit, this issue is believed to depend sensitively on the number of fermion 
species Nf. Chiral symmetry breaking is supposed to persist in the continuum limit for 
Nf < N fc, where the critical value Nf c has been estimated as anything between | and 

oc l,^ l ^l with recent estimates based on truncated Schwinger-Dyson equations yielding 
Lattice estimates have proved inconclusive due to the enormous separation of 
scales predicted by the SD approach (eg. the combination (tp^/e 2 ~ O(10~ 4 — 1CT 5 )): 

so far studies on lattices up to 80 3 suggest Nf c > lU^ 

It is intriguing that in QED3 the two classic non-perturbative phenomena are ac- 
counted for by very different theoretical approaches; confinement is attributed to instan- 
ton effects whereas chiral symmetry breaking is accounted for by self-consistent solution 
of the SD equations. An interesting issue raised by these considerations is whether com- 
pact and non-compact models have the same continuum limit, which seems probable if 
indeed monopoles are irrelevant for IR physics. Are any other continuum scenarios pos- 
sible beyond a confining gapped theory and a deconfined chirally symmetric one? How 
in fact is "confinement" characterised in the presence of dynamical electric charges? 
How does everything depend on the parameter Nf? These questions may have impor- 
tant potential a pplicati ons in condensed matter physicspEESIand have stimulated recent 
numerical work.^HZI 



2 The Lattice Model 

In this paper we will present the results of numerical simulations of a variant of compact 
lattice QED3 in which an extra four-fermion contact term has been added to the action, 
which in terms of real- valued link potentials 9^ x and iVy/2-component staggered fermion 
fields ip x , Tp x reads 

s = J2MD xy + M xy )ip y + ^J24 + P E (i-cose^) (l) 

xy x x,[i<u 



where D xy = \ E M (e ie ^ 5 yjX+( i - e ^Sy^-p), M xy = (M + § E <x ,x> °x)5xy, 
and Q^ux = A+6 UX — A+6^ x . Here 7]^ x are Kawamoto-Smit phases, a a real-valued 
scalar auxiliary field defined on the dual lattice sites x, and < x, x > denotes the set of 8 
dual sites neighbouring x. Gaussian integration over a yields an attractive four-fermion 
interaction of the form —G(t^tp) 2 , with the coupling G oc 1/ yj(3 s . Pure compact lattice 
QED3 is recovered in the limit S — > 00. 

The four-fermi term has been introduced in studies of Ad models — to enable sim- 
ulations in the massless limit M. = 0. Although in the limit f3 — > 00 it induces chiral 
symmetry breaking in its own right for some /3 S < j3 sc , in Ad the extra interaction is 
irrelevant, and should leave the continuum limit unaltered as j3 s is made large. In 3d the 
situation is less clear: (3 SC 0.25Nj defines a UV-stable renormalisation group fixed 

pointlP^ so that while we may suppose that as @ s //3 is increased any non-perturbative 
behaviour such as chiral symmetry breaking may be attributed to gauge dynamics such 
as monopoles, it remains to be checked that this behaviour is not simply associated with 
the fixed point of the 3d Gross-Neveu model. Note also that the four-fermi term reduces 
the chiral symmetry group from U(l) to Z2. 

Magnetic monopoles in the lattice model Q are identified a la DeGrand-ToussaintP^ 
for every plaquette the Dirac string content is identified via 

where G (— vr, n] and s is integer. Gauge invariant integer magnetic charges on the 
dual sites are then given by 



in? 



tfiuX^SvXx. (3) 



Since on a 3-torus the total magnetic charge J2 X ™x = 0> we it useful to define 
the average magnetic charge magnitude per site m = V~ x J2 X l^sl as a measure of 
monopole activity. 

We have performed simulations with Nf = 4, M. = 0, using a hybrid molecular 

dynamics algorithm. Fig-^sets the scene. We show both S = (a) = (5^ 1 (ipip) , and (m) 
as functions of (3, for two four-fermi couplings S = 2.0 and f3 s = 4.0. Broadly speaking 
we see that chiral symmetry breaking is approximately restored at large 0, the transi- 
tion occuring at a coupling (3 c {f3 s ) which grows stronger as (3 S is increased. There are 
corresponding peaks in the scalar susceptibility at the same locations (see Fig.|5]below), 
suggesting that a true phase transition occurs. There is not yet enough data on different 
volumes to make a definitive statement about the nature or order of this transition; since 

it has proved so difficult to identify Nj c from simulations on finite systemsP^lwe should 
not exclude the possibility of a first order transition to a high-/? phase where chiral sym- 
metry is still very slightly broken. What we can say, though, is that the fermions in the 
weakly coupled phase > j3 c are much lighter. We also see that (m) decreases monoton- 
ically with [5, with both ('ipip} and (m) becoming "small" at /3 w 1.5 for j3 s = 2.0. This 
correlation between chiral condensate and magnetic charge density was first observed by 

Fiebig and Woloshyn^l and more recently by Fiore et a$r^ The correlation looks less 
convincing, however, for S = 4.0, since changing (3 S has relatively little effect on (m). 
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Figure 1 : Scalar condensate E and magnetic charge density (m) vs. /3 for two values of j3 s 



3 A Fresh Approach 

As Adriano once memorably reminded us, ascribing properties such as confinement to 
a system because of the presence of monopoles is a little like claiming a metal is super- 
conducting simply because it contains electrons. A more sophisticated characterisation 
of the vacuum is needed. Unfortunately, here it is impossible to construct a disorder pa- 
rameter signalling a broken symmetry associated with monopole condensation, as in Ad 
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abelian lattice gauge theory, because as mentioned above in 3d monopoles are not par- 
ticles. Instead, we will follow an approach introduced in studies of 3d non-abelian gauge 

theoriesP^J an d look at the difference in free energy and its derivatives, and in particular 
their dependence on lattice volume, when a magnetic test charge is introduced. Physically 
this is akin to probing the properties of the resulting magnetic field in the presence of both 
dynamical monopoles and dynamical electrons: is it screened, indicating persistence of 
confinement, or can its effects operate over arbitrarily large distances? 

First let us review the construction, beginning with the quenched case. The free energy 
of the system is implicitly given by 

Z = e-? F = J V9^e- S , (4) 

with as usual 

S = (3 (1-cosB^). (5) 

Assuming that the action in the magnetic charge q ^ sector has the same form, we can 
rewrite the same expressions for a system in the presence of an external monopole: 



(6) 



Sm = P ^2 i 1 ~ cos ®i*»x) ■ ( 7 ) 
with 0^,, defined below. The difference in free energy we are looking for is then formally 

AF = F M -F=hog^-. (8) 

P 



Because of the difficulty in measuring partition functions in lattice simulations, we 



defim 



,24 



/' = ^(/^f) = (,S'.w),s, f -(.S').s, «9) 



The subscript indicates the action with respect to which the average is taken. To deter- 
mine (Sm}s m m independent simulation using the action Sm instead of S is needed, 
implying a duplication of computational effort. By integrating p we could in principle 
determine AF. However, unless we are interested in its precise value, the general fea- 
tures of the physics can be extracted directly from p: eg. we might expect it to display a 
sharp peak in the critical region if monopoles have something to do with any phase tran- 
sition. We can also monitor other observables, such as Am = (m)s M — {Tn)s, 
etc. 

For the form of the action Q, a naive proposal is to use the expression 

®tiux = ®(j,vx ~\~ -B^ ux , p,, v = 1, 2, 3 , (10) 

where -B Mra = A+b„ x — A+6 MX is the plaquette built with the link field 

- q fxn 

b(r) = -— ^— , (11) 

2 r(r — r ■ n) 

which is the vector potential of a Dirac monopole located at r = carrying q units of 
magnetic charge with string along h. The action in this form would have nice properties, 
like the correct continuum limit and invariance under gauge transformations on the field 

b. However, if we adopt this form, we cannot use periodic boundary conditions. In fact, 

it is easy to see that a simple link shift 9^ x i— ► 8^ x + b^ x will reabsorb b into the partition 
function, so Zm = Z and no information on AF can be extracted. 

Our proposal is to keep standard boundary conditions for all dynamical fields, but to 
use free boundary conditions on the magnetic field. In other words, the lattice topology 

seen by b is not that of an hypertorus, but that of a standard cube. This ensures a net 
magnetic flux exiting the box, and hence by Gauss's law a non-zero magnetic charge. 
This "mixed topology" requires extra care for the plaquettes on the bottom face (i.e. 
with a coordinate equal to 0) of the lattice: since we have to take account also of the 
contribution of the external field on the top face, these plaquettes enter the action twice, 
with different B^ v . The action then reads 

Sm/P = zZi 1 -^ cos (€V + B ll) ~ \ cos(6 Mi , + B d ° v wn )\ (12) 
+ (1 - cos (G M „ + B^)) , 
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Figure 2: Quenched results: (left) In p and (right) ln(Am/p) vs. /3 for various lattice sizes. 



where P is the generic plaquette, B is the boundary of the lattice and the superscripts up 
and down refer respectively to the top and bottom faces of the lattice. 

When fermions are introduced, the interaction with the the external field entails the 
following replacements in D xy : 

giOfix > ^{O^x+b^x) _ ej-®^ x Q^ x (13) 

With B' the boundary without the edges, £ the edges and M. the bulk without the bound- 
ary, and denoting a generic link by I = (x, p) we have: 

A link in the edge must be identified with three other parallel links, here called upl, up2 
and upYl. Note that the ( couplings on the boundary are not unitary. 



(14) 



4 The Quenched Theory 

We first illustrate these ideas in the quenched model, where we "know" the answer. The 
left panel of Fig. |2 presents In p(/3) for lattices ranging from 16 3 up to 32 3 , and shows 
a crossover from strong to weak coupling behaviour at j3 « 1.5, whose location is vol- 
ume independent. The behaviour at weak coupling can be accounted for by Polyakov's 
model of the QED3 vacuum as a dilute neutral monopole plasma. The screened Coulomb 
potential is 

-Mr 

V(r) oc (15) 

r 

withal 

M 2 = 4vr 2 /3exp(-/3(l7) = 4vr 2 /3£ (16) 



where £ is the monopole fugacity. We can use this effective potential to estimate the 
excess action in the presence of a test magnetic charge. The Lagrangian density C oc 
H 2 (r), with H oc -W. Therefore 



d 3 rC(r) oc 4tt 



dre 



-2Mr 



M + 



1 



1 2 



(17) 



The UV contribution dominates the integral; since Ma ~ exp(— Cj3) we deduce that at 
weak coupling 

p „ lV 2Ma oc exp(-i7T^Pe- c ^ 2 ), 
a 



(18) 



and thus that p depends only weakly on (3 in the continuum limit. The data of Fig. |2] 
support this, with the proviso that even on 32 3 there is no sign of the saturation as L — » oo 
predicted by (fTTl . 

With a few more assumptions it is possible to estimate Am. The density of m (m) 
is proportional to £e ±gV ( r \ where g = 2ir^/f3 is given by Dirac quantisation. With the 
(surely too crude) approximation that m and m densities are uncorrelated we can write 



Am oc £ / d 3 r[cosh(2-jry/pV(r)) - 1]. 



For large r the integrand may be approximated as 2it 2 PV 2 , so 



Am oc ^pe-°^ 2 exp(-4ir^f3e- c ^ 2 ) 



Am 
P 



oc y/fie 



-CP/2 



(19) 



(20) 



The right panel of Fig. |2] shows that the ratio Am / p predicted in d20b is qualitatively 
correct for and moreover volume independent as L — * oo, suggesting that the 

simulation is able to probe the dilute plasma characteristic of the continuum limit. Note 
that the /3-axis has been reparametrisedf' 



2 In 



-i 



(21) 



to match the Wilson to the Villain form of the lattice QED3 action, for which the Polyakov 
picture should be more accurate. Unfortunately quantitative agreement is less good; our 

data yield a value for C smaller than the analytic prediction^^ and the positive curva- 
ture at large j3 is hard to explain, probably requiring a better treatment ofm-m correla- 
tions. 



5 Results 

We now present results obtained from O(10 6 ) HMD simulation trajectories of each of the 
actions <Q) and (112113b on 16 3 , 20 3 and 24 3 lattices, using both /3 S = 2.0 and /3 S = 4.0. 

The left panel of Fig. |3] shows p((3) for the three different volumes at S = 2.0. 
The curves peak at (3 1.3, with no significant evidence for the position of the peak 
changing with volume. As in the quenched case, p increases with volume and there is 
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Figure 3: p vs. /3 for /3 a = 2.0 on various lattices (left), and for different p s on 24 s (right) 
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Figure 4: Am vs. /3 for various f3 s on 24^ (left), and \i vs - /9 f° r 9 — 0, 1 on various lattices (right) 



no sign of saturation; however, there is also no sign of a divergence, either at a particular 
j3* , signalling a phase transition, or for a range of values of (3, signalling an infinite free 
energy for the q ^ sector and hence a phase in which magnetic charge is confined. In 
the right hand panel we compare the 24^ results of the two fermionic models with the 
corresponding quenched data. There is remarkably little difference between j3 s = 2.0 
and j3 s = 4.0; the main effect seems to be a shift of the peak towards stronger coupling 
compared with quenched, a generic result of electric charge screening due to virtual 
fermion pairs. The peak location is independent of j3 s , unlike the apparent location of 
the chiral transition f3 c seen in Fig.^ So far there is no reason to interpret the p(/3) data 
as describing anything other than a crossover from strong to weak coupling, as in the 
quenched case. 

The left hand panel of Fig. |4] tells a similar story for Am. This quantity probes the 
cloud of virtual mm pairs produced by the plasma to screen the external magnetic charge. 
Once again there is a well-defined peak, this time at a slightly stronger coupling, and still 
well to the left of the quenched peak. An interesting feature is that the two fermionic 
curves only differ significantly in shape in the range 1.0~/3Sl.4 - Fig. ^ confirms that 
in this region the fermions are massive for the j3 s = 2.0 theory but massless (or at least 
much lighter) for j3 s = 4.0. This suggests that light fermions are able to suppress the 



production of virtual mm pairs to some extent. There is, however, no evidence that the 
large-/? behaviour differs from the quenched theory. 

Finally, the right hand panel of Fig.|4]shows the chiral susceptibility xi = 9 (ipip) / dM. 
(or at least its disconnected component) versus f3, separately in both "normal" q = and 
"monopole" q = 1 vacua on the three different volumes. The peaks are consistent with 
the chiral phase transition at (3 C ~ 1.4 seen in Fig. Q but it is interesting that chiral 
symmetry restoration seems to occur at a slightly weaker coupling in the presence of an 
external monopole. This can be understood semi-classically; the energy of a spin-singlet 
ipip pair is lowered in the presence of a magnetic field since the constituents have parallel 
magnetic moments, so pair condensation is promoted. The vacuum therefore contains 
more ipip pairs in the vicinity of the external monopole, and the resulting inhomogeneity 
in (ipip) increases xi- Fig.|4]shows, however, that at least at large (3 there is no significant 
increase of this effect with volume, suggesting that the spatial region over which pair 
enhancement takes place is finite in extent, implying in turn that the external monopole's 
field is still Debye-screened by the mm plasma even in the presence of light dynamical 
fermions. 

In summary, the data so far are consistent with the scenario that the monopole plasma 
survives the introduction of dynamical fermions, even once chiral symmetry is appar- 
ently restored. We have identified a region of crossover from strong to weak-coupling 
behaviour similar in nature to that seen in the quenched theory, and there is no sign of 
any phase transition or singular free energy in the magnetic sector. Moreover there is no 
obvious relation between the chiral transition and the crossover. At this stage, therefore, 
we favour the scenario of Ref. 10 over that of Ref. 9, with all the usual caveats about the 
need for a more quantitative data analysis, better data in the critical regions, and better 
understanding of the volume scaling. . . 

While we have at least a provisional answer to the question of monopole survival, the 
nature of the continuum limit of compact QED3 remains intriguing. Is chiral symmetry 
truly restored, at least for Nf > 4? If so, we have an example of a theory which is both 
"confining" and chirally symmetric! Perhaps this can be understood in the sense of a 
finite-ranged interaction between conserved fermion currents; we can speculate whether 

the continuum limit coincides with those of the xU<p3 fermion-gauge-scalar modelP^or 

even the 3d Thirring modeJ^at their UV fixed points. An alternative scenari is that 
chiral symmetry remains broken in the continuum limit, but at a scale far too small to 
be detected on currently available lattices. In this case the distinction between compact 
and non-compact theories persists as the issue of whether the photon is gapped or not. 
What is clear is that QED3, apparently the simplest of gauge theory models, continues to 
tantalise us even after more than 15 years of numerical study. 
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